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I. INTRODUCTION 

Two objects that are ubiquitous in the Quantum Optics lit- 
erature are the qubit (that is, a quantum two-level system) 
and the quantum harmonic oscillator (oscillator for brevity). 
From a mathematical point of view they are two of the sim- 
plest quantum systems, and yet they have been employed to 
model accurately a wide variety of physical phenomena. As 
a notable example, composite qubit-oscillator systems have 
been very successful tools for modelling and understanding 
the interaction between matter and quantized radiation 1 1 1 . 

Recently, however, substantial theoretical and experimen- 
tal efforts have been directed in a somewhat opposite direc- 
tion. Rather than using these two objects merely as theoretical 
tools to model "naturally available" systems, current research 
focuses on engineering quantum systems such that their be- 
haviour resembles that of qubits or oscillators as closely as 
possible. In addition, it is desirable to have accurate exper- 
imental control over these systems and their mutual interac- 
tions. This shift of perspective is partly due to the emergence 
of research subjects related to Quantum Technologies, such 
as Quantum Information |2|, and partly due to interest in the 
foundations of Quantum Mechanics. 

In this spirit, experimental platforms such as nanomechan- 
ical resonators |3|, ion traps |4|, cavity QED |5| and circuit 
QED |6| have been capable of demonstrating a variety of 
qubit-oscillator models in the lab. Not only can these systems 
reproduce fundamental light-matter interactions on a larger 
scale and/or in completely different settings, they are also able 
to explore these interactions in parameter regimes that are un- 
available at the microscopic level [7 |, and even simulate qubit- 
oscillator Hamiltonians that have a qualitatively different form 
as compared to the "natural" interactions. This increased flex- 
ibility suggests the possibility of studying Hamiltonians that 
are alternative to e.g. the standard Jaynes-Cummings model 
of Quantum Optics |8 |, with a focus on their applicability to 
Quantum Technologies and to foundational studies of Quan- 
tum Mechanics. 

Here we are concerned with one such model. We consider a 
particular form of interaction Hamiltonian between qubits and 



quantum harmonic oscillators, whose closed system dynam- 
ics generates qubit-controlled displacements. That is, the time 
evolution resulting from our considered Hamiltonian shifts the 
average values of the position and momentum of the oscilla- 
tor, by an amount that depends on the internal state of the 
qubit. This simple mechanism can establish strong quantum 
correlations between the two parts of the system |9 |, which 
in turn can be exploited in many applications that are of gen- 
eral interest for both technological and fundamental reasons. 
These include the preparation of non-classical oscillator states 
IfTOll . the preparation of entangled states that are robust against 
thermal noise lfTn[T2ll . the reconstruction of the quantum state 
of single oscillators fT3^/T4| and oscillator networks 1 15 1, the 
realization of quantum communication tasks such as quantum 
teleportation or remote state preparation |9|, the violation of 
Bell inequalities 1 16 1, the thermometry of quantum oscillators 
ifTTl . the measurement of geometric phases |18|, the prob- 
ing of oscillator coherences |19| and many others. Even a 
proposal for full quantum computation has been put forward 
which relies on the same mechanism EOl . In view of these 
and other proposals, the need has arisen to identify suitable 
experimental set-ups where the model can be implemented, 
and to investigate quantitatively how decoherence affects the 
realization of the desired tasks. 

In this paper we give a contribution towards answering 
such questions. We show how to realize the proposed qubit- 
oscillator model in several realistic experimental platforms, 
namely nanomechanical systems, ion traps, cavity QED and 
circuit QED, and we solve the master equation resulting from 
the inclusion of common decoherence mechanisms affecting 
both qubits and oscillators. This can be of general usefulness 
to anyone interested in the study of the proposed system. We 
then provide two relevant examples in which our techniques 
are applied to study the influence of the thermal environment 
on physical quantities of interest. First we study the time evo- 
lution of qubit-oscillator entanglement, arguably the key re- 
source in many of the applications described above. As a sec- 
ond example, we investigate the capability of the system to 
produce oscillator states with negative Wigner functions. 

The paper is organized as follows. In Section [ll] we intro- 
duce the model of interest, in particular presenting the gen- 
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eral forms of the Hamiltonian and master equation that will 
be considered throughout the paper In Section [III] we provide 
a review of some experimental platforms where the model can 
be implemented, and for each considered system we include 
a detailed discussion of the necessary steps to obtain a mas- 
ter equation in the desired form. In Section |IV] we show how 
the master equation can be solved in its general form, and in 
particular we derive analytical expressions which are valid in 
many cases of interest. In Section [V] our techniques are ap- 
plied to study the dynamics of qubit-oscillator entanglement 
and the preparation of non-classical oscillator states, while in 
Section rvTl we draw our conclusions. 



Comparing the above expression with the definition of the dis- 
placement operator D(a) = go-a'-o-'a^ ^^jj i-gcast the time 
evolutor in the compact form 

U(t) = D((TMt)), (7) 

a{t) = -/ r dsgis)e"\ (8) 
Jo 

The operator of Eq. (j7]i can be interpreted as a qubit-controlled 
displacement operator, that is, a displacement operator for the 
oscillator, where the sign of the parameter [±a'(f)] is deter- 
mined by the eigenvalue of the qubit operator 0-3. 



II. THE MODEL 



B. Open system dynamics 



We consider a qubit with ground state \g) and excited state 
|e), and an oscillator with annihilation operator a. The two 
interact via a Hamiltonian that conserves the excitations of 
the qubit. We suppose that such interaction is bilinear with 
respect to the qubit operator 0-3 [where (0-1,0-2, 0-3) is a right- 
handed tern of Pauli matrices] and the oscillator position, and 
we allow the coupling strength g{t) to be a generic function of 
time. The Hamiltonian thus reads 



(1) 



where v is a frequency that does not need to coincide with the 
natural frequency of the oscillator 



A. Closed system dynamics 

For the time being, let us assume that the system is closed, 
i.e., that its interaction with the external environment is negli- 
gible on the timescales of interest. Under this assumption, the 
density matrix obeys the Liouville equation li2TII 



p(f) = -i[H(t),p{t)l 
This equation is formally solved by 

Pit) = U{t)p{Q)U{t)\ 



(2) 



(3) 



where p(0) is the initial state of the system, while U{t) is 
the unitary time-evolution operator, which verifies the time- 
dependent Schrodinger equation: 



U{t) = -iH{t)U(t), 
t/(0) = 1. 



(4) 
(5) 



Solving Eq. Q is in general difficult, since it involves the cal- 
culation of a time-ordered exponential. In our case, however, 
the time ordering operation just yields a global phase factor 
which can be ignored for our purposes ifTSl . The time evolu- 
tion operator can thus be expressed as 



U(t) = exp 



-iO-3 



J dsg(s) (c 



(6) 



So far we have neglected losses in the system, which al- 
lowed us to obtain a unitary time evolution. A realistic model, 
however, has to take into account the weak but inevitable cou- 
pling of the system to its surrounding environment. Here, we 
shall consider a widely applicable model of Markovian deco- 
herence for both the qubit and the oscillator |22 |. To take into 
account environmental effects, we write down a master equa- 
tion by modifying Eq. (|2| as follows: 



-i[H(t),p]+Lp + Qp. 



(9) 



The effect of decoherence is included through the following 
Lindblad terms, which are responsible for the non-unitarity of 
the dynamics: 

L = "^{Na + ima] + '^NaDla^l (10) 

Q = jiN, + l)D[cr-] + ^NgD[cr^] + jDlcr^l (11) 

where £)[A], A being a generic operator, yields the Lindblad 
form 



D[A]p = 2ApA^ - AUp - pA^A. 



(12) 



In Eqs. (10 1 and (111, the parameters k, Fi and F2 quantify 
the eff'ective coupling strength of the system to the environ- 
ment. Typically, k and Fi are associated to the exchange of 
excitations between system and environment, while F2 to ex- 
tra dephasing mechanisms affecting the qubit. Note that we 
have assumed that the environment acts independently on the 
two constituents of our system, and in particular, that the de- 
coherence terms in the master equation (|9]l are not affected by 
the variations in the qubit-oscillator coupling. This treatment 
is consistent provided that the density of states of the environ- 
ment is approximately constant in the frequency range rele- 
vant to the system [see e.g. Ref. [23 1 for a similar issue arising 
in driven solid state qubits]. The environment is assumed to 
be at thermal equilibrium at temperature T, hence A^^ and A^^ 
indicate the number of thermal excitations associated to the 
natural frequency of the oscillator and qubit respectively (de- 
pendence on the temperature T has been left implicit). When 
the environment is modelled as a collection of harmonic oscil- 
lators, its number of thermal excitations at a given frequency 
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is given by the Bose-Einstein distribution A^^ - (e'^J^'^ - 1)"' 
where ajj, for j = q, a, is the natural frequency of either the 
qubit or the oscillator. 



III. EXPERIMENTAL IMPLEMENTATIONS 

Before attacking Eq. (j9]l, let us review a few experimental 
setups that can be used to implement our model. In the fol- 
lowing examples, our aim will be to identify physical systems 
where the Hamiltonian ([T]i and the master equation (|9| can be 
reproduced with good approximation. Depending on the par- 
ticular setup under examination, it can be the case that even 
though the "bare" Hamiltonian is not in the desired form, it 
is still possible to obtain a Hamiltonian equivalent to Eq. ([T]l 
through appropriate unitary transformations and approxima- 
tions. When this is the case, the applicability of the deco- 
herence model of Eq. (|9]l has to be considered carefully, as 
one has to take into account also the effect that the mentioned 
transformations have on the decohere nee terms of Eqs. (lOi 
and ( [TT| i. We shall see that, for the cases considered in this 
paper, it is always possible to obtain a master equation in the 
same form as in Eq. (|9|, although the meaning of the various 
decoherence parameters might be different as compared to the 



standard scenario provided in Section II B 



A. Flux qubit coupled to a nanomechanical resonator 

It is known that our considered model can be implemented 
by coupling nanomechanical resonators to superconducting 
qubits Il24ll . In particular, we shall focus here on the system 
put forward by Fei Xue ef. al. |25|. In their theoretical pro- 
posal, the coupling between a superconducting flux qubit ll26l 
and the fundamental flexural mode of a nanomechanical os- 
cillator |3 1 is controlled by the amplitude of a transverse mag- 
netic field, resulting in a Hamiltonian 

s 6 .1. / t\ 

^^flux(f) = zO"3 + Tzo-i + uja' a + g(t)o-i [a + a' j , (13) 



git) = T]Bit), 



(14) 



where s is the level splitting of the superconducting qubit, 6 
the tunnelling energy between the two levels, ai the mechan- 
ical oscillation frequency, rj a constant that depends on the 
specific system realization, and B{t) the externally controlled 
magnetic field amplitude along an appropriate direction. A 
remarkable feature of this setup is that the three parameters 
s, 6, git) can be tuned independently to each other, so that the 
Hamiltonian of Eq. ([T]l is easily obtained by operating in the 
regime 6 ~ 0, and considering an interaction picture trasfor- 
mation with respect to the free term |cr3 + cja'a. As a typi- 
cal set of parameters |25 -28J, we can take s - In x lOGHz, 
oj ^ 2kx lOOMHz, T] ^ 2nx 0.8MHz/mT. From the forni 



of Eq. ( 14 1, we see that both the magnitude and sign of the 
coupling can be controlled over time by tuning the external 
magnetic field. Assuming a maximum magnetic field inten- 
sity Bmax - lOmT, we get that the magnitude of the coupling 




FIG. 1 : Sketch of a flux qubit coupled to a nanomechanical resonator, 
as described in Ref. |25 1. In this picture, the flux qubit is given by 
a superconducting loop containing three Josephson Junctions (indi- 
cated with the symbol "X" in the picture), while the nanomechanical 
resonator is embedded in the loop, oscillating perpendicularly to the 
plane of the circuit. The magnetic field Sy, parallel to the plane of the 
circuit, and perpendicular to the direction of the mechanical oscilla- 
tion, controls the strength of the qubit-oscillator interaction, while 
the perpendicular field Bx can be used to tune the free Hamiltonian 
of the qubit. See Ref. L25J for more details. 



strength g(t) can vary between zero and g^ax - 2;7r x 8MHz. 
For the decoherence parameters of the system we can take 
K ^ 0.3(//s)"\ ^\ - ^2 - 2.5(yL(s)"', while the thermal ex- 
citations of the environment are given by A^^, ^ 0, Na - 20 
at a temperature T ~ lOOmK. In units of the mechanical 
frequency, this parameter range yields a maximum coupling 
strength gmax - 0.08w, a qubit decay rate Fi ^ 4 x IQ^^cl) 
and dephasing rate y ^ O.Olw, a mechanical coupling to the 
environment k ^ 5 x 10 '^w and a heating rate kNo - O.Olw. 



B. Trapped Ions 

Ion traps f4\ constitute a rather mature experimental plat- 
form, which is especially suited for the implementation of our 
model. In this case two electronic levels of the ion provide 
the qubit, while the motion of the ion inside the trap can be 
modelled as a quantum harmonic oscillator with good approx- 
imation. For simplicity, we assume that we can consider only 
the motion of the ion along one of the trap axes f4U), say x. 
To couple the internal levels of the ion to its motion along the 
X axis, we drive the system resonantly with a standing laser 
wave, in such a way that a node of the electric field coincides 
with the centre of the trap ||29l . Correspondingly, the total 
Hamiltonian of the system reads 

Hionit) = ^cr3-i-wfl'a-i-2i2(?)cri sin |^^xo (a + cos st, (15) 

where s is both the frequency of the ion electronic transition 
and that of the driving laser, ll) the ion motional frequency 
along X, Q(f) the Rabi frequency of the driving standing wave, 
with corresponding wavenumber k, xq the ground state spread 
of the ion along the x axis. We allow the Rabi frequency to 
vary over time in a controlled way, provided that this variation 
is slow compared to the transition frequency |j41J s. Taking an 
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interaction picture with respect to Hq - |cr3, we get 

H[^^{t)-u)a^a+Q.{t) (cri +cr^e^""+cr"e""^"") sin [^xo(a+fl"'")] . 

(16) 

Note that s can be in the optical or infrared domain, while for 
the mechanical frequency a typical value is w ^ 2:7r x lOMHz. 
Therefore, if we are interested in the dynamics on a timescale 
of u), we can consider a Rotating Wave Approximation (RWA), 
and neglect the fast oscillating terms at frequency ±2e |30|. In 
addition, since typical ground state spreads for trapped ions 
are of the order of xq ~ 1 lOnm, while the wavenumber 
is of the order k ~ lO^m ' for optical transitions, it follows 
that = kxo <K 1. If in addition the motional state of the 
ion is in the Lamb-Dicke regime 77, -y/((fl + a"'")^) <K 1 , we can 
approximate sin [z/y (^a + a'^j ^ rj^-ia + a' ). As a result, the 
Hamiltonian of Eq. ( jlS] ) is well approximated by 

^k,n(0 ^ « + 8(t)o-i {a + fl"'") , (17) 
git) = TjMt), (18) 

which can be transformed into Eq. (jT]) via an interaction pic- 
ture with respect to wa'a and a qubit rotation cr^ — > cr3. Note 
that there is a trade-off between the accuracy of the linear ap- 
proximation of the sine, which improves as rj^ gets smaller, 
and the strength of the qubit-oscillator coupling, which is di- 
rectly proportional to 77^. For example, considering a Lamb- 
Dicke parameter 77 ^ ^ 0.05 can guarantee linearity of the cou- 
pling within a reasonable range of motion without weakening 
too much the coupling strength g{t). By choosing a stable tran- 
sition of the ion, decoherence of the qubit can be neglected, 
since the electronic levels involved can have lifetimes up to 
several seconds. Note however that the Rabi frequency scales 
as n oc y/ir, where I is the laser intensity and F the de- 
cay rate of the excited state, so that stabler transitions require 
higher laser intensitites to be driven |4|. Rabi frequencies up 
to ilmax - 2n X IMHz have been reported for dipole forbid- 
den transitions with lifetime ~ Is |31|, which together with 
our choice rj^ ^ 0.05 would allow \g(t)\ to vary between zero 
and a few percent of the motional frequency a>. If the driving 
standing wave is realized by means of an optical cavity (oper- 
ating in the semi-classical regime), this figure can be increased 
further The time variability of the coupling can be preserved, 
provided that the cavity response time is small compared to 
the motional timescales of the ion. As anticipated, the avail- 
ability of long internal lifetimes allows us to set Q = in the 
master equation (|9|. On the other hand, motional decoher- 
ence is usually relevant, and typical heating rates in ion traps 
at room temperature are of the order kNo ~ lO"* lO'^s"'. 
Thus, at a temperature T ^ 300K, we have Na - 6 ■ 10^ 
for <x> - 2ji X lOMHz, and we may take k ^ 0.02s"'. In 
units of the mechanical frequency, we have a: ^ 6 x 10"^ w, 
kNq - 2 X IQr^o), and by considering a cavity-assisted driv- 
ing we may assume g^ax - 0.05w. Typically, the Lamb-Dicke 
regime in which our model is valid is reached by laser-cooling 
the ion to an excitation number much lower than N„. 



FIG. 2: Sketches of trapped ions. Left: several ions confined in a 
linear trap. Right: a laser standing wave is used to couple tiie motion 
of an ion along tlie x axis witli two of its internal levels. In order to 
realize tiie desired Hamiltonian tlie center of tiie ion motion has to 
coincide with a zero of the driving field \29L 



C. Cavity QED 

Cavity Quantum Electrodynamics (QED) has long been 
the paradigmatic quantum optical setup to investigate simple 
qubit-oscillator models |5 1. The qubit is provided by two elec- 
tronic levels of an atom or ion (atom for brevity), while the os- 
cillator represents a single confined mode of the electromag- 
netic field inside a high-finesse optical or microwave cavity 
(cavity mode for brevity). The atom can be either trapped 
inside the cavity, so that it interacts with the cavity mode con- 
tinuously, or it can fly through the cavity and interact with 
its confined mode only for a limited time. The interaction 
between a two-level atom and a single-mode electromagnetic 
field is described by the Rabi model 



^^Rabi(0 = |o-3 + Loa'^a + g(t)cri (a + a"'") , 



(19) 



with s the level splitting of the atom, oj the frequency of the 
cavity mode, g{t) the interaction strength between atom and 
cavity mode at time t. The latter is usually well approximated 
by considering ^(f) = ^0 while the atom is within the effective 
volume of the cavity mode, and ^(f) = otherwise. Despite 
its rather simple form, the Rabi Hamiltonian of Eq. (19 1 is 
in general difficult to treat analytically 1321 . However, if the 
atom and the cavity are resonant or near-resonant (s ^ w), and 
their coupling strength is weak compared to the uncoupled fre- 
quencies [go E, a>), it is a good approximation to describe 
their interaction via the Jaynes-Cummings Hamiltonian |8|: 



Hjc(t) = -0-3 + Loa^'a + g(t) {, 



a' cr + ao- 



(20) 



The Hamiltonian of Eq. ( 20 1 has innumerable applications. 



nevertheless it leads to a very different dynamics from what 
we are pursuing in this paper One possible way of obtain- 
ing a Hamiltonian of the form ([TJ in cavity QED is to add a 
strong external driving to the atom |10|. Here, we adopt the 
complementary strategy of strongly driving the cavity, an ap- 
proach that might be useful in situations where direct strong 
driving of the atom is not feasible. Consequently, we modify 
the Hamiltonian as 



Hdri 



-(T3 -H cija^ a + go (ao-^ + a ' cr ^ + 



+ n[ae-''^^"" + a^Q''^- 



(21) 
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where Q is the strength of the cavity driving, operated at reso- 
nance with the qubit frequency e, (f is the phase of the driving 
field at time f = 0, and we considered the atom to be inside 
the cavity during the time of interest, so that g{t) - gQ. Going 
to a rotating frame at the driving frequency s, we get 



6a^a + go (flo-+ + a"' "cr") + Q (ae"''^ + a"' e'^) , (22) 



where 5 = w - e is the detuning between atom and cavity. We 
see from this form of the Hamihonian that the "equiHbrium 
position" of the oscillator in phase-space has been shifted. To 
find the new equilibrium position, it is necessary to take into 
account decoherence. The master equation reads 



(23) 



where L and Q keep the same form of Eqs. ( 10 1 and (111 also 
in the rotating frame we are considering. We now multiply 



the master equation of Eq. ( 23 1 by the displacement operator 
D(q'o)' on the left, and D(ao) on the right, so that we reaUze 
the bosonic translation a — > a + q-q- Imposing that the trans- 
formed master equation does not contain terms proportional 
to {ap,a^ p,pa,pa^), we get the value of ckq corresponding to 
the new phase-space centre of oscillation: 



Qe'^ 



0!Q - -- 



'7 



(24) 



For simplicity, we choose the driving phase ip such that q-q > 
0. Correspondingly, the transformed master equation reads 



P = -'[^d'rive'P]+^P + 2P' 



where the new Hamiltonian is given by 



^a + gQ\acr + a' cr 



(25) 



(26) 



Note that Q' = g^ao is an effective driving strength for the 
atom. We are now very close to getting our desired Hamil- 
tonian. Taking an interaction picture with respect to the last 



term in Eq. (26 1, we obtain a modified master equation 



P^-i[H'^L(t),p]+Lp + Q'(t)p, 



(27) 



^drive(0 = + ^goo-i (fl + fl') -H ^god-2(t) (fl - a') , (28) 
Q'(t) = y (A^, + l)D [i(cri - icr2(t))] + 

^N,D[^icTi+iiT2(t))] 



(29) 



where 0-2(1) - cr2Cos2Q.'t - cr3sin2Q'f and 6-3(1) - 
0-3 cos 2Q'f + 0-2 sin2Q'f. If the strength of the effective driv- 
ing is large enough, such that ll42lll43l Q.' » go,S, we can 
apply again the RWA, and neglect the fast-oscillating terms 
at frequency ±2Q' |30|. Applying this approximation to the 
Hamiltonian is straightforward, while some care has to be 
taken in applying it to the operator Q'(t), as it contains terms 
proportional to cos^ 2Q'f or sin^ 2Q.'t, which are fast oscillat- 
ing but do not average to zero. One can verify that the correct 



approximation is then 
D[l(cri±icr2(t))] - \D[cri] + ^©[ctj] + ^-Dlcr^] (30) 
D[cr3(t)]^^-D[cr2] + ^D[cr3], (31) 

resulting in 



Q' ^^D[(rn + jn[(T2] + ^D[(r3], (32) 

n = Ti (a^, + i) , 72 = ^ [n, + '2) + r2. (33) 

One may also notice the relation 

^(D[cr2]+D[cr3])^D['^]+D[^]. (34) 

Taking a further interaction picture with respect to 6a'' a, and 
a qubit rotation ct] — > 0-3, we finally get the master equation: 



P = -iWcqed(f), P]+Lp+ Qp, (35) 

//cqed(f) = g(t)(r3 {ae-'" + a^e'") , (36) 
Q^^{D[cr^]+D[cr-])+jn[cr3]. (37) 



where g(t) - go/ 2 if the atom is inside the cavity, g(t) - oth- 
erwise. By comparing Eqs.([TT| and Eq. (37i, we can see that 
in the latter case the heating and cooling rates of the qubit are 
equal, and that the role of decay and dephasing rates has been 
swapped. Another difference from the previous examples is 
that here the tunability of the coupling strength is restricted to 
switching between zero and a fixed value, essentially by con- 
trolling how much time the atom will spend inside the cav- 
ity. On the other hand, the atom-cavity detuning 6, which in 
this case plays the role of the oscillator frequency, might be 
modified from experiment to experiment by applying a Stark 
shift to the atom[44J. By considering circular Rydberg atoms 
flying through superconducting microwave cavities Il33l . we 
can have coupling strengths g - go 1 2 ^ 2n x 25kHz, to- 
gether with low decoherence rates k ^ l(ms)"' for the cavity 
and Fi ^ 0.03(ms) ' for the atoms ||5]- Considering cavity 
modes of frequency o) - 2n x 50GHz operated at a temper- 
ature T ^ 1.3K as reported in |33|, we have A^^ ^ 0.15 for 
the atoms, resulting in y\ - 2y2 - 0.04(ms)"'. For the cavity 
mode the effective thermal occupation of the environment is 
slightly higher due to room temperature photons leaking in- 
side the cryogenic chamber: A^^ ^ 0.8. At resonance (6 - 0) 
the strong driving condition ao » 1 implies that the cavity 
driving should be stronger than its decay rate: Q » /c. As 
here the coupling g is fixed, it makes sense to re-express the 
other parameters in terms of this quantity: ^ 6 x IQr^g, 
kNo - 8 X IQ'^g, yi = 272 - 2 X 10""*^, and a total dephasing 
rate y = 71 h- 72 - 3 x 10""^^. 



D. Circuit QED 

In circuit QED f6l, the Rabi model of quantum optics is 
realized"on-chip", by coupling superconducting qubits (the 
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qubit-oscillator coupling ifTTl [T9ll and a zero-temperature en- 
vironment for the qubit lfT3l . Here we discuss how to solve 
the master equation in its most general form. 



A. C-Matrix representation 



FIG. 3: Sketch of a cavity-QED setup. Two internal levels of an atom 
or ion are coupled to a mode of the electromagnetic field, confined 
between two highly reflective mirrors. The thick arrow on the left of 
the picture represents the external driving laser, necessary to obtain 
the desired Hamiltonian. 



artificial two-level atoms) to the quantized excitations of mi- 
crowave resonators (i.e. microwave photons). In this archi- 
tecture, light-matter interactions typical of cavity QED can 
be simulated in parameter ranges that are unavailable at the 
single-atom level |7|. Regimes where the Jaynes-Cummings 



model of Eq. ( 20 1 is valid are routinely achieved in circuit 



QED |34|, so that the same methods we discussed in Sec- 
tion III C can be straightforwardly applied here to obtain the 
Hamiltonian of Eq. (36 1. Recently a novel superconducting 
qubit architecture has been proposed by the Houck group in 
Princeton, which allows for the independent tuning of both 
the qubit energy splitting and its coupling strength to the mi- 
crowave photons |35|. Tunability of the coupling strength go 
between a negligible value and 27rx45MHz has been demon- 
strated in the Jaynes-Cummings regime, with the qubit and 
the cavity at resonance, oj - s ^ Inx 5.8GHz |36|. The same 
group demonstrated that this tunability can be achieved on 
the timescale of few ns [37 J . Consequently, this architecture 
should allow the realization of the Hamiltonian of Eq. ( [36| , 
where g(t) can be an experimentally controlled function of 
time. Typical decoherence parameters are /c ^ 7.5(yus)"', 
Fi ^ 0.6(yus)"', r2 - 0.5(jus)"' while the system is operated at 
the low temperature T ^ 25mK, yielding Nq - Na ~ 0, and 
ji ^ 0.3(;us)-i,r2 ^ 0.65(ps)-i 136, 37 1. 

Even though what we have described until now applies to 
the Jaynes-Cummings regime, the versatility of the tunable 
qubit architecture should in principle allow to operate in a 
regime where the qubit energy splitting e is negligible as com- 
pared to the qubit-cavity coupling git). Then, the Rabi Hamil- 
tonian alone would be sufficient to provide the required dy- 
namics, without the need of an extra driving of the cavity. In- 
deed, by neglecting the qubit frequency in Eq. ([T9]l, we get 



^^Rabi(f) - <wfl"' fl + g(t)cri [a + a"'") , 



(38) 



which is equivalent to Eq. ([T]i, up to an interaction picture 
transformation. 



IV. SOLVING THE MASTER EQUATION 

The master equation Q has been solved in the literature 
in several special cases of interest. These include a constant 



To attack the master equation we adopt a phase- 
space representation for the oscillator degrees of freedom 
ifTTl [T3l [T5l |T9l . In particular, we shall simply extend the def- 
inition of the oscillator characteristic function |22 | to a hybrid 
qubit-oscillator system. Note that the same idea can be ap- 
plied to other phase-space representations, such as the Wigner 
function |[n][T9l. We define the C-Matrix as 



(39) 



where tr„{ ) indicates partial tracing over the oscillator Hilbert 
space. Many properties of the characteristic function are eas- 
ily extended to the C-Matrix. In particular, the knowledge of 
the C-Matrix is equivalent to the knowledge of the quantum 
state p, due to the inverse relationship p - \ ^ d^f^x(fi')D(-ff). 
The master equation (|9]) is equivalent to the following partial 
differential equation (PDE) for x, where the complex param- 
eters /3 and p* are treated as independent variables Il22l : 



X = -ig(t) {e'^'d^ - e-'^'d^.) [cr,,x]+ 



In Eq. (Boll, Q has the same form as in Eq. ( 1 1 



(40) 



is the 



anticommutator, while the operator responsible for the deco- 
herence of the oscillator has now the differential form: 



(41) 



with A = Na + h . Since Eq. (|40|l only involves the qubit op- 



erators 1X3, 0"+, o" , it becomes convenient to decompose x in 
terms of the eigenstates of 0-3, as 



X(P, t) = Xee(J^, t)\e){e\ +XgM t)\8)<8\^ 
+ XeMO\e}{g\+Xse(/^,t)\g}{e\. 



(42) 



As a result, we get the following system of PDEs for the C- 
Matrix elements: 

Xeg^2ig(t)(e-<''df^.-e''df,)Xeg+£Xeg-7Xeg, (43) 

Xge ^-2ig{t){t-"'dp, -t"'dp)Xge+£Xge-yXge, (44) 
Xee^ig{f){t-"'l}^e'"l3*)Xee+-LXee-TcXee+ThXgg, (45) 

Xgg ^^g{t){^-"'l3+e''nXgg+LXgg-'^hXgg+'^cX..- (46) 

In the above equations, y - ri(A^fy + 1/2) -i-r2 can be inter- 
preted as the total dephasing rate of the qubit, while r^, = 
^i{Nq + ^),^h - ^iNq are respectively its cooling and heat- 
ing rates. Note that for the special cases of cavity-QED 
and circuit-QED, the decoherence parameters in Eqs. ( |43| l- 
( |46l ) are instead given by 7 = yi +72 - jT] (Nq + jj + r2, 

rc=r,,=72 = ^(A^,+i)+r2. 
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B. Off -diagonal terms 



We can see that Eqs. (43 1 and (44i are in decoupled form, 
and in fact they can be solved analytically for any choice of 
the parameters. The corresponding solutions are given by 



where 



Jo 



fdi^We'^sinh^i, 
smh ff Jo 2 



Jo 



(47) 
(48) 



(49) 
(50) 
(51) 



Some important applications such as the reconstruction of the 
oscillator state ifTTl [131 only make use of the two C-Matrix 
elements Xeg,Xge^ so that in this case the above solutions al- 
low to analyse exhaustively any parameter regime where our 
model is valid. On the other hand, it is obvious that the matrix 
elements Xei>^Xi>e alone are not enough to describe the com- 
plete dynamics of the system, and in general the remaining 
elements ;^ff,e,A'i;g will be needed. 



C. In-diagonal terms 



Eqs. (45 I and (46 1 represent a system of two PDEs coupled 
to each other via the heating and cooling rates of the qubit, 
and their solution presents more difficulties. We shall show 
here how it is possible to reduce them to a system of ordinary 
differential equations (ODEs) in t, for each value of the com- 
plex parameter /3. Analytical solutions are presented for the 
special case in which the heating rate vanishes. 



1. General treatment 

We shall start to build our solutions for the C-Matrix ele- 
ments by considering the uncoupled problem obtained 
by setting T/, = in Eq. ( |45| ) and Ff. = in Eq. ( |46l ). It can be 
checked directly that this is solved by 



;e,.(/?,0 = e-^'-^('-^') 



Jo 



-''W+^«*^;r,g(y8e^',0), 
^+^(fr-aW/3^^^(^e-f/^0), 



(52) 
(53) 

(54) 



Then, we can use the uncoupled solutions to write down an 
ansatz for the coupled problem, in the form 



Xgg(J3, t) = xggd^, t)Q>gg(l3e^-', t), 

Xeeifi, t) = Xeeiji, O^eeiP^''^' , t). 



(55) 
(56) 



After some manipulations it can be verified that the above ex- 
pressions represent the exact solution to Eqs. ( |45] l and (|46|, 
provided that, for each fixed value of y6, ^ggiji, t) and <i>ee(J3, 
as functions of time solve the following system of ODEs 

*..Oe,0=r,^:^^e<r'.^^)-2^^'^'«/^'-'«-^io,,(/3,0, (57) 

Xggip,^}) 

'i'..(;8,0 = r,^^e(r'^'>-2ef'[i(o-Wl,5^^(^^,), (58) 

with initial conditions <S>gg(j3,Q) = <S>ee(J3,0) = 1. In general, 
the above system of ODEs can be solved numerically for each 
value of jS, which provides a net computational advantage as 
compared to attacking directly the system of PDEs (|45| and 
( |46l ). We stress that, once the funcions '^e,g(fi, t) are available, 
the substitution fi — > jSe^i' has to be performed before these 
are used to build the full solutions to Eqs. (|55]l and ([56)1. 



2. Perturhative solution 

In many cases one is interested in interaction times such 
that F/,?, Fc? «: 1. Then, a reasonable approximation to the 
exact solution to Eqs. ( |57] i and ( |58] l can be given at first order 
in time-dependent perturbation theory: 



%gm - 1 + F, 



Xeelfifi) 
XggiPfi) 

xgg(m 



Jo 



^^^2.2'lM.s>p--Aism^ (59) 



(^ee(/3,t) ^ 1 + r/, ^^"'";;,'^; r dse2'=''[-^<^>*'^'^<^'^*l. (60) 
Jo 



Xee(J3,0) 



Vanishing heating rate 



In the special case F^ ^ the equations decouple, and in 
particular Eq. ( |58| ) simply gives (tgj, - 0, hence Off(/S, f) = 1. 
Then, Eq. ( [ST] ) is straightforwardly integrated as 



a),,Oe,f) = 1 + F 



Xee(J3,0) 

^XggiPfi) 



Jo 



(61) 



In applications where F/, is interpreted as a physical heating 



rate [such as Sections III A and IIIB| , this corresponds to con- 
sidering an environment for the qubit that does not possess 
thermal excitations: N^j ^ 0. Note that while assuming this, 
we will still allow the possibility of having a thermally excited 
envkonment for the oscillator: A^,, 0. This can be the case 
when the frequency of the qubit is several orders of magnitude 
higher than that of the oscillator (this is the case e.g. in flux 
qubits coupled to nanomechanical oscillators or ion traps, as 
discussed above). 



V. APPLICATIONS 

Having shown how to derive our model in realistic experi- 
mental setups and discussed how to solve its dynamics in the 
most general form, we move on to apply our techniques to two 
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practical examples that can be of general interest: the study of 
entanglement in mixed qubit-oscillator states and the prepara- 
tion of states with non-classical Wigner functions. To keep the 
mathematical complexity to a bare minimum, we shall con- 
sider a simple instance of our model: we take the Hamiltonian 
of Eq. ([T|l with a constant coupling profile g(t) = go and a 
vanishing oscillator frequency v = 0. This can be the case. 



e.g., in cavity QED and circuit QED setups (Sections IIIC 
and |IIID| i. In order to avoid the technicalities associated with 
the qubit heating and cooling rates (see Section [TV C| l we set 
Tft = Fj, - 0, and consider a qubit decoherence of pure de- 
phasing: y (in fact, when considering the data presented 



in Sections III C and III D one might neglect qubit decoher- 
ence altogether, as a first approximation). 



A. Dynamics of qubit-oscillator entanglement 

It Ref. 19] it was shown that the mechanism of qubit- 
controUed displacements is able to maximally entangle the 
qubit-oscillator system, even when the oscillator is initially 
in a highly mixed state. Our aim is to investigate how this 
prediction is affected by the presence of the environment. We 
assume that the system is initially in the separable state 



po = |+X+|®Pth, 



(62) 



where |±) - (\e) + \g)) / y/l are the eigenstates of cri, and 
Pth - Zin/'nl«)(n| is the state of the oscillator when it is at 
thermal equilibrium with its environment, so that |n) are the 
Fock states and p„ = N"/{Na + l)"'^^- We now consider the 
time evolution of such a state according to the master equation 



(|9]l, and the techniques of Section IV can be used to provide 
the solution for any f > 0. The initial conditions for the C- 
Matrix elements that have to be used are Xj(P,0) = ^e"'^^'", 
for j - e,g, +, -, and we recall that A = Na + l/2. After some 
manipulations it can be shown that the C-Matrix elements at 
time t are given by 



;^,,06,f)=ie-^^-2.oP-H. 



\2 



(63) 
(64) 
(65) 
(66) 
(67) 

(68) 



where for brevity we have left implicit the functional depen- 
dence of the various parameters. These correspond to the den- 
sity matrix 



p = i [|e)<e|D(ao)Pth£>"'"(«o) + \g){g\Dhao)p,uD{aQ)\ 
^ [|e)<g|D(ao)Pth£'(ff()) + \g){e\DHaQ)p,^D\ao)] , (69) 

as it can be checked directly by substituting the above expres- 
sion in Eq. ([39|). In Ref. |9|, the entanglement of the state p 



was studied in the limit of large displacements ([ffol » VA). 
Here, on the other hand, we cannot in general rely on this ap- 
proximation, and we have to attack the problem without fur- 
ther assumptions on the parameters entering the state p. We 
shall quantify the qubit-oscillator entanglement via the nega- 
tivity ll38]|||45l . We start by calculating the partial transpose of 
p with respect to the qubit: 

pT' = i [|e)<e|£)(ao)Pth£>^(ao) + \g){g\D\aQ)pi,D{aQ)\ 

[|g)(e|D(ao)Pth£>(«o) + k)te|£)"f(ao)Pth£)^(ao)l • (70) 

If we find a vector such that {\Ij\p'^i\\Ij) < 0, then the state 
p must be entangled f38l. To this end, following f9l, to each 
Fock state \m) we associate the test state 

\ilf,n) = -^^{\e)D\ao)\m)-\g)D{ao)\m)). (71) 

Defining q,„ - (i/'/,,|p'^'|i/'m), and using the fact that the test 
states are orthonormal {{ij/mVl'm') - S,nm'), we can find a lower 
bound for the negativity of the form 



Nip) > S^ip) = 2 ^ \q„ 

?,„<o 



(72) 



Note that the quantity is in itself a valid entanglement 
witness. In particular, a positive value signals the presence of 
entanglement, while a value of 1 implies that p is a maximally 
entangled qubit-oscillator state. Evaluating S^(p) is straight- 
forward once the expectation values q,„ entering Eq. ( |72| i are 
known. These can be formally expressed as 



q,,. 



(73) 



where L„, is the Laguerre polynomial of order m. The typi- 
cal behaviour of S^ip) is shown in Fig. |4|a), where we can 
see that the presence of qubit-oscillator entanglement can be 
witnessed in a range of temperatures and interaction times. 



B. Preparing non-classical superposition states 

It is known that a Hamiltonian of the form ([T]l can be used 
to prepare non-classical superposition states of the oscillator 
(see e.g. ifTOl [13 24|). Indeed, it is easy to check that by 
preparing a state |+) (gi |0), applying the time evolution oper- 
ator of Eq. (|7|, and measuring the qubit in the basis |±), the 
oscillator is left in a superposition of coherent states of the 
form \a-t) cc \a} + \ -a}. The non-classicality of these states is 
typically associated to their Wigner function assuming nega- 
tive values ||39l . Thanks to our techniques, here we can inves- 
tigate how realistic decoherence and thermalization mecha- 
nisms influence the generation of such states, and in particular 
the negativity of their Wigner functions. Before proceeding, 
we recall that given an oscillator state pa with characteristic 
function ;^f a (yS) = tra{Di/3)pa], its Wigner function is given by 
d 



(74) 



9 




FIG. 4: (a): Lower bound for the negativity of the state p, quantified 
via the entanglement witness 23^. (b) Plot of "W, indicating the ca- 
pability of the system to produce states of the form p_, with negative 
Wigner functions. All quantities are plotted as functions of the oscil- 
lator average thermal excitation A'^ and rescaled interaction time ggt, 
the decoherence parameters being set to k = y = O.Ol^o- The two 
quantities in plots (a) and (b) show a similar qualitative behaviour, 
reinforcing the intuition that qubit-oscillator entanglement is a key 
resource for engineering non-classicality in the oscillator. 

Our starting point is again the separable state ( |62] i, which is 
the natural generalization of \+) igi |0) to the case where the 
oscillator environment is at nonzero temperature. Following 
the usual strategy, we let the system evolve for a time t, thus 
obtaining the qubit-oscillator state p, as given in Eq. ( |69l ), and 
finally we measure the qubit in the basis |±) so that, depending 
on the outcome, the oscillator is projected onto the states 

p^^(p^r\±\p\±}. (75) 

In the above equation, Pj- - tia{{+\p\+}} is a normalization 
constant, corresponding to the probability of obtaining the 
outcome |±). For simplicity we shall restrict to the outcome 
|-), as the Wigner function of the state p_ exhibits the maxi- 
mum negativity (if any) in zero. Correspondingly, we define 
the following quantity 

^ = 7rP_max{0,-W(0)), (76) 

where W(a) is the Wigner function of p_. Note that 'W en- 
codes information about the non-classicality of the state p_, 
as well as the probability of preparing it by measuring the re- 
source state p [Eq. (|69])]. It may thus be interpreted as an 
indicator of the system's capability to produce non-classical 
oscillator states. The normalization constant in Eq. (|76]l is 
chosen such that 'W varies between and 1 for the considered 
class of states (see below). In particular, a value of 1 signals 



the optimal situation in which a state with W{0) — -l/n can 
be produced with probability P_ = 1 /2. Explicit calculation of 
the quantities involved in Eq. ( |76l ) yields 

1 _ g-4A|Q'ol"-n' 

P = ^ , (77) 

e^-e-" 
W(0) = . (78) 

A typical plot of 'W is shown in Fig. |4|b), where we may 
notice a similar qualitative behavior as compared to the en- 
tanglement witness [Fig. |4|a)]. Indeed, in this context it 
is generally perceived that the possibility of generating an os- 
cillator state with negative Wigner function is directly linked 
to the presence of qubit-oscillator entanglement in the state p 
of Eq. ( |69l ). The quantitative analysis presented here seems 
compatible with this intuition. 



VI. CONCLUSIONS 

We have studied a qubit-oscillator model in which a Hamil- 
tonian generating qubit-controlled displacements competes 
with a thermal Markovian environment. We have shown how 
the model can be experimentally meaningful, in particular 
for nanomechanical systems, ion trap, cavity QED and cir- 
cuit QED. We have shown how the dynamics of the system 
can be solved with minimal numerical effort, and analytically 
in many relevant special cases. Finally, we have applied our 
model and techniques to show that the system can become 
entangled and show non-classicality of the oscillator Wigner 
function even in the presence of rather strong thermal noise. 
Our techniques can be useful to study the experimental fea- 
sibility of many applications that rely on the mechanism of 
qubit-controlled displacements. 
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